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Determination of S17 from 7 Be(d, n) 8 B reaction: 
CDCC analyses based on three-body model 
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The astrophysical factor Sn for 7 Be(p, 7) 8 B reaction is reliably extracted from the transfer re- 
action 7 Be(d, n) 8 B at E — 7.5 MeV with the asymptotic normalization coefficient method. The 
transfer reaction is accurately analyzed with CDCC based on the three-body model. This analy- 
sis is free from uncertainties of the optical potentials having been crucial in the previous DWBA 
analyses. 

PACS numbers: 24.10.Eq, 25.60.Je, 26.65.+t, 27.20.+n 



The solar neutrino problem is a central subject in the 
neutrino physics |l]] . The major source of the high-energy 
neutrinos observed by solar neutrino detectors is 8 B pro- 
duced by the 7 Be(p, 7) 8 B reaction. The astrophysical 
factor 5i7 for the reaction, however, is one of the most 
poorly determined reaction rates in the standard solar 
model; the latest recommendation for the factor 5i7(0) 
at zero energy, based on recent direct measurements, is 
19±2 eV b [|, i.e., 10-20 % error exists. This is far 
from our goa of determining 6*17(0) within 5 % error 
required in order to determine the neutrino oscillation 
parameters: the mass difference between v e and and 
their mixing parameter. The main difficulty in the direct 
measurement comes from ambiguities of determining the 
effective target thickness of the radioactive Be beam. 
Thus, indirect measurements of 5*17(0) are expected to 
be essential for determining Si 7(0) accurately. 

The transfer reaction 7 Be(<i, n) 8 B at low energies is an 
example of such indirect measurements; once the asymp- 
totic normalization coefficient (ANC) of the overlap func- 
tion between the 7 Be and 8 B ground states is determined 
from the data of the transfer reaction, Si7(0) can be 
accurately derived from the ANC, as long as the reac- 
tion is peripheral 0. The reaction has been measured 
at E = 7.5 MeV and the ANC is extracted with the 
distorted wave Born approximation (DWBA) [Q. The 
Si 7 (0) obtained with the ANC is 27.4 ± 4.4 eV b, lead- 
ing to inconsistency with the recommended value. Al- 
though the reaction is found with DWBA to be indeed 
peripheral ^ |(| , distorting potentials used in DWBA are 
quite ambiguous j^, [|, which prevent one from extract- 
ing Si7(0) accurately. In particular, uncertainties of the 
d- Be optical potential bring about large errors for S17, 
typically 30 % in magnitude. The origin of the large am- 
biguity of distorting potentials is that these are derived 
from proton and deuteron optical potentials for different 
targets and/or energies. 



In the present paper, we analyze 7 Be(<i, n) 8 B at 7.5 
MeV with the three-body model, p + n+ 7 Be, assuming 
7 Be to be an inert core. This treatment of the system 
is justified by Refs. [Q, |^], where core excitations of 7 Be 
are shown to be negligible. An advantage of this analysis 
is that we do not need the ambiguous c?+ 7 Be optical po- 
tential in the entrance channel and the ambiguous n+ 8 B 
one in the exit channel as shown below. 

The three-body dynamics in the entrance channel 
are explicitly treated by means of continuum-discretized 
coupled-channels (CDCC) method || |Io[|, the theoreti- 
cal foundation of which is given in Ref. [Jll|. This theory 
has been established as a method of solving the three- 
body system with good accuracy, and extensively applied 
for various reactions || [l2|]. Previous CDCC calculation 
showed that explicit treatment of breakup channels is es- 
sential in describing deuteron induced reactions Q. The 
CDCC thus provides a precise description of the wave 
function in the entrance channel, i.e., c?- 7 Be system. 

The effective Hamiltonian for the entrance channel 
based on the three-body model contains the optical po- 
tential between N (p or n) and 7 Be. Data of the neutron 
elastic scattering are available for target 7 Li, the mir- 
ror nucleus of 7 Be, at 4 MeV, approximately half the 
deuteron energy considered here |13| . First the n- 7 Li 
potential is determined accurately from the data. The 
potential is then used as an input in CDCC calcula- 
tion for deuteron elastic scattering on 7 Li at 8 MeV, and 
the numerical result is compared with the experimental 
data |Q. This is a good test for the neutron optical 
potential determined above, which is only an input in 
CDCC calculation for the entrance channel. As for the 
exit channel, the three-body dynamics are treated with 
the adiabatic approximation, after testing its accuracy 
with CDCC calculation. The ANC is then obtained with 
reliable distorted wave functions in both the entrance and 
exit channels. 

The transition amplitude for the transfer reaction, 
based on the three-body model (p + n+ 7 ~Be), is 



'Electronic address: kazu2scp@mbox.nc.kyushu-u.ac.jp 



Tr, 



(-) 



( + ) 



> 



(1) 



2 



The three-body wave function is a solution of the 

Schrodinger equation (H — E)^\ + ' = for the three- 
body Hamiltonian 



Hi = K ap + K dBc + V Bp (r np ) + U, 



(2) 



with U = J7 n Bc(i - nBc) + ^pBc(i - pBo)- Here r XY is the coor- 
dinate of nucleus X relative to nucleus Y. The potential 
V np is the interaction between n and p, U p b c (U n Bc) is 
the proton (neutron) optical potential for the target 7 Be 
at half the deuteron incident energy, and K np and K d Be 
show kinetic energy operators for two-body systems de- 
noted by the subscripts. The nuclear part of U p b c is 
assumed to be the same as U n Bc- The Coulomb part 
of t/pBe is treated approximately by replacing the co- 
ordinate 7pBe by rdBc- Effects of Coulomb breakup of 
deuteron, not included in the usual treatment above, are 
found to be quite small in the present system. 

The wave function \E'j + ' > is obtained with CDCC, that 
is, by solving the Schrodinger equation in a model space. 
In CDCC, deuteron breakup states are classified with 
linear and angular momenta, k and £, and truncated into 
< k < fc max and < £ < £ max , respectively. The k 
continuum [0, k max ] is further divided into bins with a 
common width A. The total wave function is expanded, 
in terms of the deuteron ground state and the discretized 
breakup states, into 



(+) _ 



~<Mr n p)Xi( r dBc 



(3) 



where 0o is the deuteron ground state and <pi is the i- 
th discretized breakup state obtained by averaging con- 
tinuous breakup states in the i-th bin. The coefficient 
Xi represents a center-of-mass motion of n-p pair in 
the i-th state. Inserting this form into the three-body 
Schrodinger equation leads to a set of coupled differen- 
tial equations, 



(Ei - K dBe ) Xi =J2 Fi 



ij Xj 



(4) 



with Fij =< <f>i\U\4>j > rn and Ei — E — ei, where 
is an intrinsic energy of the i-th n-p state. The coupled 
equations are soluble, since they have a compact kernel 
in its integral equation form. The precise formulation of 
CDCC is shown in Ref. [^o| . The present model space is 
fcmax = 1-7 fm- 1 , t = 0,2 and A = 1.7/40 fin" 1 . The 
CDCC solution converges at these values, as the model 
space is enlarged. 

The exit channel is also treated in the three-body 
model, that is, the exit channel wave function ' is 
determined by the three-body Hamiltonian 

Hi = K pBc + K nB + VpBc(r P Bc) + U nBc (r nBc )- (5) 

In the three-body model, 8 B is treated by the two-body 
(p+ 7 Be) model with the potential VpBc- The spectro- 
scopic factor S cxp in Eq. ([!]) is introduced by taking 



account of the incompleteness of the model. It should be 
noted that Hf does not contain V^p since the interaction 
is already treated as a transition operator in Eq. ([!]). In 
general, the distorting potential U nBc between an outgo- 
ing neutron and Be differs from the corresponding one 
in Hi, since an outgoing neutron has a different velocity 
from an incoming deuteron in the (d, n) reaction. 

In principle, Hf allows transitions between the ground 
and continuum states of the p+ 7 Be system in the exit 
channel. However, effects of the transitions may be sim- 
ply estimated with the adiabatic approximation, since 
the ground state of 8 B has a binding energy (0.137 MeV) 
considerably smaller than an energy (4.18 MeV] of out- 
going neutron. Following Johnson and Soper |lq] , we re- 
place the Hamiltonian of the p+ 7 Be system by the bind- 
ing energy of 8 B. Errors of the adiabatic approximation 
are estimated with CDCC in the exit channel scattering, 
8 B(n,n) 8 B, at E = 4.18 MeV. The breakup effect itself 
is not so large and errors of the approximation are less 
than 3 % at forward angles. The Johnson-Soper approx- 
imation leads to a simple form ^[ ' = XnB&pRc, where 
</> P Bc is the wave function of 8 B in its ground state and 
X^b i s * ne w ave function of outgoing neutron distorted 
by the potential J7 n Be(' , nB "8/7), where the zero-range ap- 
proximation is made to the transition amplitude (fil) . We 
discuss the use of this approximation below. It should be 
noted that in DWBA x[\b is determined by the elastic 
scattering of neutron from 8 B, of which no measurement 
has been done so far. On the other hand, the three-body 
model approach can avoid this difficulty, as in the en- 
trance channel. 

The transfer reaction 7 Be(d, n) 8 B is calculated with 
the zero-range approximation with its finite-range cor- 
rection pH[ ], The integration over r p B G in T is made up 
to a large value 40 fm, since the transferred proton is very 
weakly bound in 8 B. The finite-range correction for the 
transition amplitude (|l|) including deuteron ground and 
breakup channels is straightforward; the resultant cor- 
rection for the i-th channel keeps the standard form by 
regarding Fu and as the potential and the intrinsic en- 
ergy of the entrance channel. This prescription is tested 
by doing finite-range DWBA calculation for the deuteron 
ground channel which is a main component of the transi- 
tion amplitude ([!]). The result of the prescription above 
agrees with that of the full finite-range calculation within 
2 % error at forword angles. 

The wave function 4> p b c is calculated with four types 
of VpBo 17, [l8| |l^, pCp . We determine the spectroscopic 
factor S'cxp comparing the calculated 7 Be(<i, n) 8 B cross 
section with the experimental one, for every type of the 
four potentials. The astrophysical factor 6*17(0) at zero 
energy is then obtained from the S'cxp with 



Si 7 (0) 



1 h 2 
0.026 



(6) 



in the ANC method ||, where b is defined with the 
Whittaker function W as 4> p bo{i~ p Bc) = frW^pBe) at r P B e 



3 



larger than the range of the nuclear force VpBe between 
p and 7 Be. It should be noted that if the reaction is pe- 
ripheral, ANC should be stable against the change of the 
8 B internal wave functions, i.e., that of the parameter set 
of single-particle potentials. Thus, one can estimate the 
error of the ANC calculation from the deviation of ANC 
(or Sn) with four different models of 8 B above. 
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FIG. 1: Results of the optical potential search for neutron 
elastic scattering at 4.08 MeV from 7 Li. Experimental data 
are taken from Ref. Jig) ], 
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7 Li(n,n) at 4.26 MeV 
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FIG. 3: Results of the optical potential search for neutron 
elastic scattering at 4.26 MeV from 7 Li. Experimental data 
are taken from Ref. pi)] . 
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FIG. 4: The calculated cross sections for 7 Be(d, n) 8 B at 7.5 
MeV with Soxp = 1.0 (solid line) and 0.849 (dashed line), 
compared with the experimental data M. 
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FIG. 2: Comparison between the CDCC calculation and the 
experimental data Q for 7 Li(d,d) 7 Li at 8.0 MeV. 



Figure 1 shows the result of the optical potential search 
for neutron scattering at 4.08 MeV from 7 Li, the mirror 
nucleus of 7 Be. The _resultant potential shows a good 

The optical potential is then 



agreement with data 1 1 3 



applied for deuteron scattering at 8.0 MeV from 7 Li. The 
CDCC calculation with the potential again gives a good 
agreement with data (lj] at angles 9 < 60°, as shown in 
Fig. 2. The potential is shown to be reliable especially 
at the forward angles; we use it in Hi and obtain the 
proper wave function for the entrance channel. As for 
the exit channel, on the other hand, we need the optical 
potential of the n+ 7 Li scattering at E = 4.18 MeV to 
determine Hf. The data is available at E — 4.26 MeV 
that is closest to the proper energy. Figure 3 shows the 
result of the optical potential search for the scattering. 
The resulting potential well reproduces the data {Hf at 



TABLE I: Parameters for the optical potentials between neu- 
tron and 7 Li at E n — 4.08 (a) and 4.26 MeV (c) corresponding 
to the initial and final channels for 7 Be(d, n) 8 B at Ed~ 7.5 
MeV, respectively. The sin gle particle potential between p 
and 7 Be in 8 B of Kim et al. (b) is also shown. 
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ro 


ao 


w d 


Ti 




Vso 




Oiso 


a 


46.57 


2.07 


0.49 


0.82 


1.87 


0.22 


5.50 


1.15 


0.50 


b 


32.12 


1.54 


0.52 








8.24 


1.54 


0.52 


c 


60.97 


1.47 


0.58 


0.31 


3.57 


0.22 


9.0 


2.39 


0.55 



Figure 4 shows our result for 7 Be(d, n) 8 B cross section 
using V^,Be of Kim et al Jl7|., compared with the exper- 
imental data JiJ. The solid and dashed lines represent 
the calculated results with S c ^ p = 1.0 and 0.849, respec- 
tively. Parameter sets of the optical potentials at 4.08 
and 4.26 MeV and the single-particle potential of Kim et 
al. used in the calculation are listed in Table 1 together. 
At forward angles 6 < 60°, the calculated cross section 
well reproduces the data H] with the spectroscopic fac- 



tor S exp = 0.849, leading to Sn(0) = 21.36 eV b. It was 
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found that deuteron breakup states play important roles 
not only in determining distorting potentials but also in 
the transfer process. In fact, when the deuteron breakup 
components is set to zero in , the resultant transfer 
cross section is reduced by 10 % at the forward angles 
and, more seriously, the angular distribution cannot be 
reproduced correctly. The components are obviously not 
included in the framework of the standard DWBA. We 
thus conclude that the three-body model approach is in- 
evitable. 



TABLE II: Results of SWp, b, and Sn(0) with different 8 B 
single particle models; see Eq. (0) for definition. 







Scxp 


6 


Sn(0) 


Kim et al 




0.849 


0.809 


21.36 


Tombrello 


m 


0.882 


0.784 


20.87 


Robertson 


fit 


0.864 


0.794 


20.93 


Esbensen and Bertsch |2C|] 


1.097 


0.700 


20.67 



We show in Table 2 the list of calculated S cxp , b, and 
5i7(0) for different VpBe- One sees that the calculated 
values of 517(0) are almost consistent, which shows that, 
as mentioned above, the present reaction is peripheral 
and the ANC method works well. Taking account of 
the theoretical errors of ANC, the adiabatic approxima- 
tion (AD) and the finite-range correction (FRC) and of 
the systematic error of the experimental data [|4| on the 
7 Be(d,n) 8 B cross section, we obtain S 17 (0) = 20.96±g'| 
(ANC) ±0.63 (AD) ±0.42 (FRC) ±2.7 (expt) eV b, in 
consistent with both the recommended values Q and the 
recent result 5*17(0) = 22.3 ± 1.2 eV b of acccurate direct 



measurement [ p2| . 

In summary, the present analyses based on the three- 
body model are free from uncertainties of the optical po- 
tentials in both the entrance and exit channels which 
were the most essential problem in the previous DWBA 
analyses. The deuteron breakup process in the incident 
channel is significant for determining 517(0) within 5 % 
error required from the neutrino physics. We then con- 
clude that the three-body model approach is essential 
and necessary. The present analyses provide a precise 
value 5 17 (0) = 20.96;j;3(theor) ± 2.7(expt) eV b; the 
theoretical ambiguity of £17(0) is 6-7 % slightly beyond 
the required accuracy. However, we can reduce the theo- 
retical error to ~ 2% coming from ANC only, if we do the 
finite-range calculation with accurate \&j and ' de- 
rived by CDCC. Such full-fledged calculations are highly 
expected. In the present analyses, however, the experi- 
mental error (13 %) is even larger than the theoretical 
one. It is expected that the peripheral properties, es- 
sential for the ANC method, become insufficient, as the 
incident energy increases [pj. Thus, accurate measure- 
ments on 7 Be(d, n) 8 B and Li(d, <i) 7 Li at about a few 
tens of MeV and on 7 Li(7V, 7V) 7 Li at the half the corre- 
sponding deuteron incident energy are highly expected; 
the proton and deuteron elastic scattering are necessary 
to determine the nucleon optical potential accurately. 
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